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Intracavity and external third order correlations in the damped nondegenerate parametric oscil- 
lator are calculated for quantum mechanics and stochastic electrodynamics (SED), a semiclassical 
theory. The two theories yield greatly different results, with the correlations of quantum mechanics 
being cubic in the system's nonlinear coupling constant and those of SED being linear in the same 
constant. In particular, differences between the two theories are present in at least a mesoscopic 
regime. They also exist when realistic damping is included. Such differences illustrate distinctions 
' between quantum mechanics and a hidden variable theory for continuous variables. 
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I. INTRODUCTION 



Local hidden variable theories have been extensively compared to quantum mechanics over the last seventy or so 
years Most comparisons between the two have investigated whether or not quantum mechanics is equivalent 

to a local hidden variable theory. Much evidence indicates that it is not. Many results in quantum mechanics have 
been found that are incompatible with all local hidden variable theories [§-0. Most of these results have involved 
^C) ' idealized undamped systems. However, all experimental systems encounter damping. Thus, it is interesting (and 
more realistic) to compare quantum mechanics and local hidden variable theories in damped systems . This paper 

■ compares quantum mechanics and one local hidden variable theory (SED) in such a system. 

One of the earliest works comparing local hidden variable theories to quantum mechanics was Bell's theorem [|J. 

■ It demonstrates that quantum mechanics is incompatible with all local hidden variable theories at a statistical level. 
O |' It does so by deriving an upper bound on a function of two particle correlations for all local hidden variable theories, 
, 1 1 , which quantum mechanics exceeds. Extensions of it have been formulated for large angular momentum and particle 

number systems ||||. These extensions demonstrate nonclassical behavior in a regime usually regarded as being 
purely classical. Greenberger, Home and Zeilinger (GHZ) jij have also extended Bell's work, differentiating quantum 
mechanics from all local hidden variable theories for single, as opposed to ensemble, measurements. The three particle 
GHZ theorem has an "all or nothing" quality and distinguishes between local hidden variable theories and quantum 
mechanics in a single experimental run, once three basic correlations are established. 

Comparisons between quantum mechanics and local hidden variable theories have also been made using continuous 
variables (which are discretized in formulating the comparison), such as quadrature phase amplitudes and there 
is currently much interest in this area. For quadrature phase amplitude measurements, these comparisons can have 
detector efficiencies of in excess of 99% [g. They also tend to relate more strongly to Einstein, Podolsky and Rosen's 
original EPR paradox |j| than earlier discrete variable ones. Indeed, the EPR paradox has been experimentally 
demonstrated using quadrature phase amplitudes Additionally, quantum teleportation has been achieved using 
quadrature phase amplitudes [jlO| further demonstrating the utility of continuous variables. 

One commonly used local hidden variable theory is stochastic electrodynamics (SED) |Il]JlS^| . Some authors have 
pro posed it as an alternative to quantum mechanics |H],[b5]]. Furthermore, a semiclassical approach equivalent to it 
[Q is also commonly used in parametric oscillator calculations SED consists of adding Gaussian white noise 

to classical electrodynamics. It is equivalent to truncating third order derivative terms in the quantum mechanical 
Moyal equation, a commonly used approximation [fl6| . Such terms are often negligible and thus SED reproduces many 
results of quantum mechanics ]l4| , |l5 |. However, it cannot violate Bell inequalities for quadrature phase amplitude 
measurements, and is thus distinct from quantum mechanics 0. Various authors have explicitly shown differences 
between SED and quantum mechanics ]T7j-|l9|] . In particular, it has been shown that the two theories predict dif- 
ferent transient third order correlations for the undamped nondegenerate parametric oscillator Jl7[ . It has also been 
shown that they predict different macroscopic quadrature phase amplitude correlations in the damped nondegenerate 
parametric oscillator in the steady state ftj| . 



1 



In general, differences between quantum mechanics and local hidden variable theories are reduced or eliminated by 
damping ppf . Furthermore, damping is a significant element of many realistic systems. It is thus important to consider 
its effects on differences between quantum mechanics and local hidden variable theories such as SED. However, all 
but a few of the comparisons between the quantum mechanics and local hidden variable theories referenced above 
have involved undamped systems. They are thus idealized in this respect. In contrast, damping is included in the 
calculatons in this paper. It is included to consider a theoretical model which is as realistic as possible and also to 
determine the sensitivity of differences between quantum mechanics and SED to its presence. 

This paper extends a previous comparison between quantum mechanics and SED in the nondegenerate paramet- 
ric oscillator jl?]]. In particular, it contrasts both intracavity and external moments of the two theories' in the 
same system with damping included. Expressions from both theories are compared for the intracavity moment 
( AXj (t) A(l 2 (t) A J 3 (t)) , where AXj(r) = Xj(r) — (X^t)), for i = 1,2,3, Xj(r) is quadrature phase amplitude, 
the subscripts represent different radiation modes and r is a scaled time variable. A comparison is also made for an 
analogous external moment. Both analytic iterative and numerical techniques are used to calculate moments. The 
results produced by these techniques show that the intracavity and external moments differ greatly between the two 
theories. In particular, the analytic method shows that the moments of quantum mechanics are cubic in the system's 
nonlinear coupling constant to leading order whilst those of SED are linear. The two theories are compared over 
a range of nonlinear coupling constant, damping and average initial pump photon number values. The results of 
these comparisons show a number of qualitative trends. Most importantly, quantum mechanics and SED differ in 
the situations considered with the largest particle number and damping to nonlinear coupling ratios, although the 
differences are reduced in relative size. 

Stochastic techniques are used to obtain results both for quantum mechanics and SED. The positive-P coherent 
state representation pf[ | is used to calculate quantum mechanical predictions. It is particularly well suited to the 
calculation of quantum dynamics in damped quantum optical systems when nonclassical behavior is present. It is 
able to handle arbitrarily large photon numbers. It converges quickly (in the sense of sampling error) when systems' 
dimensionless nonlinearities are relatively small, as is the case with nonlinear optical experiments. By contrast, the 
method used for SED calculations corresponds to commonly used approaches in quantum optics, where the field is 
treated as a semiclassical object surrounded by (classical) vacuum fluctuations. Both methods are used to generate 
analytic predictions and are also numerically simulated. 



II. QUANTUM MECHANICS 



This paper considers an idealized nondegenerate parametric oscillator, resonant at three frequencies u>\ and L02 
(signal and idler frequencies) and W3 = u\ + u>2 (pump frequency). It contains a nonlinear medium that couples the 
modes and converts higher energy pump photons into lower energy signal and idler ones. The system's interaction 
Hamiltonian, including linear losses, is given by 

3 

H = ihG(a\a\az — aia 2 a\) + ^ f jS] + f Ja*, (2.1) 

4=1 

where a\ andfij are creation and annihilation operators for oscillator modes, T\ andl\ are environment mode operators 
and G is a nonlinear interaction strength constant. Initially, the system has a coherent state in the pump mode, and 
vacuum states in the signal and idler modes. 

A number of quasiprobability representations exist to describe quantum states, the most famous being the Glauber- 
Sudarshan representation p2|j . It is produced by decomposing quantum density operators using a diagonal coherent 
state basis. Thus, 

p = J da 2 P(a,a*)\a)(a\, (2.2) 

where p is a density operator and P(a 1 a*) is the Glauber-Sudarshan representation. The Glauber-Sudarshan rep- 
resentation can be negative and is hence not a strict probability density function. A more recent representation is 
the positive-P representation pf| which is an actual probability density function over an off-diagonal coherent state 
basis. It further differs from the Glauber-Sudarshan representation by using a phase space of doubled dimension. 
The positive-P variables {ai,a^}, where i is a positive integer, are analogous to complex field amplitudes, with on 
and af describing a particular radiation mode. However, {a{\ and {of} are independent and hence en ^ (af)*, 
though their averages are complex conjugate and thus (a,) = (af)*. Variable averages are equal to normally ordered 
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quantum averages once the substitutions a, — > a, and — ► a] are made. For example, (aia^) = (a\ai) p , where 
(0) p denotes Tr{pO), as usual in quantum mechanics. 

Stochastic equations of motion for positive-P variables for the damped nondegenerate parametric oscillator are, in 
terms of r (time scaled by T, a typical damping constant with units of inverse time), 



— — = -7!Q!i + ga%a 3 + ^fgaiii 



da i 
~d7 
da i 

~dV 

da 2 

~dT 

da£ 
~dV 
da 3 
dr 

dr 



7 2 a 2 + gafa 3 + y/ga^^ 2 



-72<4 + guiat 



"73<4 - gula+. 



Here £1 , £2 j £i~ and £, 2 are complex Gaussian white noises with the following correlations 

(&(ri)£j(r 2 )) = S 3 -i,jS(Ti - t 2 ) 
(C + (n)£+(T 2 ))^ 3 -^(n-T 2 ) 
(^(n)ei(r 2 ))=o, 



(2.3) 



(2.4) 



where i,j = 1,2. In Eq. (2.3), ji — Tj/r, where is a damping constant for mode i with units of inverse time, 
g = G/T and r = Tt . It is assumed that G, I\ and T are real. Initial conditions are «i(0) = 0, 02(0) = 0anda3(0) = e. 
It is noted that Eq. (2.3) is only valid when boundary terms in phase space can be neglected. These are asymptotically 
small i n th e limit of short times or large damping ratios f23j| . 

Eq. 2.3 is solved using an analytic iterative method. This method treats da mpin g terms exactly, and noise and 
nonlinear terms iteratively. It involves, firstly, rewriting the equations forming Eq. ( |2.3| ) as on — — 7jQ!i+/j({aj, <X^}, t) 

or af = —Jiaf + / 4 + ({aj, a^ 1 "}, t), where i,j = 1,2,3. Successively higher order approximations for {a^r), af(r)} 
are then found using increasingly better approximations for fi and Thus, (m + l)*' 1 order terms are given by 



(rn+l), s (0)/ x , 

a\ >{t) = a\ >(t) + 



dr, exp[ 7i (n - r)]M{a^\ a+ M }, n ) 



■ fl) W = 4 { °\r) + T^dnexpb^ -r)]/+({4 m) ,a+ (m) },r 1 ), 

'ri=0 



(2.5) 



where al (t) — ak{r — 0) exp(— jk T ) and a k + (r) = ai (r)* where k = 1, 2, 3. For example, 



4 m+1) (r) = a[°\r) + / dr x expMn - r)]( fl a+ (m) (n)4 m) (r 1 ) + V^fa&fa)) 

Jti=0 



and first order approximations are 

« 4 (1) (r) 
C 4 1) ( T ) = eexp(-7 3 r 



dn exp [7^(73 - r)]^/geexp( — )&(ti) 

Ti=0 ^ 



dn exp[7i(Ti - T)] v / g^exp(-^y-)4 + (ri) 



"J (1) ( r ) = e* ex P(-73T), 



(2.6) 



(2.7) 



where i = 1, 2. 
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III. STOCHASTIC DIAGRAMS 



The iterative method of the previous section can be used, in conjunction with stochastic diagrams, to readily produce 
analytic a ppr oximations for the intracavity moments of quantum mechanics considered in this paper. Stochastic 
diagrams \BM are schematic representations of the combinatoric parts of an iterative process. They clearly lay out 
all terms produced by different orders of iteration. Fundamental stochastic diagrams appear as one of three classes. 
Those associated with initial conditions appear as straight lines, those with noise terms as straight lines with a cross 
at their end and those with nonlinear terms as straight lines containing a fork, as shown in Fig. |l|(a)-(c). Higher 
order iterative terms are represented by stochastic diagrams using either combinations of the three basic classes. For 

( 2) 

example, one of the iterative terms in a\ (r) is 



r Ti T r T2 Ti 7272 4. 

/ driexp[7i(ri - T)]gal '(n) / dr 2 exp[7 2 (r 2 - ri)]Vge^exp( ^-)S~( r 2) 



dri exp[7i(n - r)}ga { 3 >(n) / 

'ti=0 Jt 2 =0 

It combines all three basic classes and is represented by the stochastic diagram in Fig. |l|(d). All iterative terms can 
be represented by stochastic diagrams. 

Stochastic diagrams can also be used to determine the orders of iterative terms. In particular, they can be used to 
determine the orders of such terms in the system's nonlinear coupling constant g. This paper focuses on the order of 
terms in this constant. For quantum mechanics, initial value iterative terms are O(g ), noise iterative terms 0(g 1 ^ 2 ) 
and nonlinear iterative terms 0{g). Hence, lines in stochastic diagrams count as order zero, crosses as order 1/2 and 
vertices as order 1. A term's order is simply found by considering its stochastic diagram and adding a half to its order 
for every cross and one for every vertex. For example, the term represented in Fig. 0(d) has one vertex and one cross 
and thus is 0{g 2 ). A notation that denotes the order in g of a term by a superscript [n] is used in this section. 

Stochastic diagrams are now used to determine the intracavity moments of quantum mechanics considered in this 
paper. Consider all eight moments of the form (AAi(t)AA 2 (t)AA3(t)), where A_4;(t) = Ai(r) — (Ai(r)) and 
Ai(r) is either hi or aj. These are equal to the positive-P variable mom ents which replace &i and aj by a^(r) and 
af(r) respectively. Now, consider the equations that constitute Eq. ( |2.3| ). Their forms do not change when they are 
expressed in terms of — oh(t), —af(r), a 3 (r) and a^{r), where i = 1,2. From this, it follows that (Ai(r)) = (— Ai(r)) 
and hence (aj(r)) = (af(r)) = 0, where again i — 1,2. Thus, {AA\{t) AA 2 {t) AA^{t)) can be simplified to 
(A 1 (t)A 2 (t)AA s (t)). _ ' 

An approximate expression for {A\{t)A 2 {t)AAz{t)) is now obtained using the iterative method in Section || 
(and stochastic diagrams). This method can be used to produce power series expressions in g for the positive-P 
variables. These expressions can then be used to generate power series expressions in g for the moments of the form 
(Ai(t)A 2 (t)AA3(t)}. As g <C 1 in realistic systems, these power series expressions can be approximated by their 
lowest order nonzero terms. 

Figs [2] (a) and (b) show the stochastic diagrams required to determine the moments of the form 
(Ai(t)A 2 (t)AA3(t)). Naively, it might be thought that the lowest order nonzero terms from Ai(r), ^4 2 (r)and AAs(t) 
simply need to be multiplied together and the average of the subsequent product determined to calculate the lowest 
order nonzero term in (A\{t) A 2 {t) A A^(t)) . This is not always true. Sometimes, A\(t), A 2 (t) and AAs(t) are not 
necessarily zero and yet (^4i(t)^42(''")A^43(t)) is zero. For example, the lowest order nonzero terms for the positive-P 
variables in (ai(T)a 2 (T)Aa 3 (r)) are 

a f /2] ( T ) = / rfr i exp[7 < (r < -r)] A /5eexp( — ^)&(t 4 ), (3.1) 
Jn=Q 1 

where i = 1 , 2 and 

A ^t [2] ( T )= dr 3 exp[7 3 (T 3 - r)]g / dr A exp[7i(r 4 - T 3 )]^/gFexp( (3.2) 

dr 5 exp[7 2 (r 5 - T 3 )]V5^cxp( (r 5 ) (3.3) 




dr 5 exp[7 2 (r 5 - T 3 )]Vge* exp( ^-)C 2 "( T 5))- 

However, the average of their product is zero as 
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(a^\r)a^\r)Aat [2 \r)) (3.4) 



f T\ —T PT2—T (• T3 — T PT4 — T3 

„3 „* 



g ee* / / / / / drsdr^drs 

Jr x =0 Jt 2 =0 Jt 3 =0 Jt 4 =0 Jt 5 =0 

exp [71(7-1 - r)j exp( — ) exp[7 2 (r 2 - r)j exp( — ) exp[7 3 (r 3 - r)J exp[7i(r 4 - r 3 )j exp( — , 

exp[7 5 (r 5 -T3)]exp(-^)«6(T 1 )6(T2)^(T4)^(r 5 ))-(6(Ti)6(T2)) (tfMtffo))) 



(6(ri)6(r 2 )^ + (r 4 )C 2 + (r 5 )) = (6(^1)6^2)) (tffo^fo)). (3.5) 



and 



Thus, the two noise terms cancel each other and the right hand side of Eq. (|3.5|) is zero. Taking such a consideration 
into account, the moments of the form (Ai(t)A2(t)AA3(t)) are determined by carefully considering the lowest order 
nonzero terms of their constituent positive-P variables and then finding the average of these variables' products. 
Consider Fig. (p|)(a), which contains the lowest order stochastic diagrams for q^t) and af(r), where i — 1,2. The 

first diagram in it represents the initial value terms a\ ^ (r) and oq + (r) , which are zero and do not contribute to 
any moments. The second represents terms containing £i,£i~,£2 or which are not necessarily zero and thus may 
contribute to moments. Fig. 0(b) contains the lowest order stochastic diagrams for a 3 (r) and a 3 (r). In it, all terms 
represented by stochastic diagrams containing initial value lines are zero except for the O(g ) ones. This is so as these 

terms contain either oq°'(r) or a^J" (t), where i = 1,2, which are both zero. In addition, all 0(g°) terms represented 
by stochastic diagrams in Fig. 0(b) are canceled out by other 0(g°) terms. This occurs because A„4 3 (r) appears 
in the moments considered. Its two components, *4 3 (r) and (^4 3 (r)), contain the same O(g ) term and hence their 
O(g ) terms cancel each other. It follows that the only remaining stochastic diagram in Fig. |(b), which represents 
the 0(g 2 ) term containing two noise components, denotes the lowest order term in A_4 3 (r) that is not necessarily 
zero. 

The lowest order nonzero terms determined above are now used to calculate (q;i(t)q; 2 (t)Aq! 3 (t)). The lowest order 
contribution to Aa 3 (r) that is not necessarily zero Aa 3 i owest (r) is 

fT 3 =T pT4=T3 



Aa 3 lowest (t) = / dr 3 exp[7 3 (r 3 - r)]g / dr 4 exp[7i(r 4 - T 3 )]y/gl exp( — )£i(t 4 ) (3.6) 

Jt 3 =0 Jt 4 =0 ^ 

r'" =T3 72T5 

/ dr 5 exp[7 2 (r 5 - Ta,)]*/ge exp{ — )f 2 (T 5 ) 

Jt 5 =0 1 

rrs—T rT4.=T 3 7lT 4 

-( / dr 3 exp[7 3 (r 3 - r)]g I (fr 4 exp[7i(r 4 - r 3 )] v /geexp( r-)£i(r 4 ) 

/ dr 5 exp[7 2 (r 5 - T 3 )]y/geexp( — )6( T 5))- 

Jt*=o 1 



When 71 = 7 2 = 7 3 = 7, the average of the product of Aa 3 \ owes t(T) and the lowest order nonzero terms in ai(r) and 
a 2 (r) is approximately equal to (ai(T)a 2 (r)Aa 3 (r)) when g -C 1 and thus 

/ t \ r u t w e 2 g 3 exp[-3jr] r exp[ 7 T] exppyr]-, , . 

(ai(r)a 2 r Aa 3 r ) ~ = x 2t . (3.7) 

7 7 

As daggered positive-P variables are complex conjugate to undaggered ones on average, (a^"(r)a 2 f (r)Aa 3 (r)) = 
(ai(r)a 2 (r)Aa 3 (r))*. The other six moments of the form {A\{t)A2{t)AA 3 {t)) are zero to 0(g 3 ). As does 
(ai (t)« 2 (t) Aa^ (r)), they all have two 0(g 3 ) terms that cancel each other. To explain such behaviour in general, 
the following argument is given. These other six moments can be rewritten as 

{Ai{t)A2{t)Mt)) ~ {Ai(t)A 2 (t))(Mt)), (3-8) 

where it is understood that the moments in which A\(t) — ai(r), Ax{t) = a 2 (r) , A3 (r) = a 3 (r) and 
Ai(t) = a+ (t), „4 2 (r) = (t) , A3 (r) = a^(r) are excluded. All 0(g 3 ) terms in the six moments of the form 
(*4i(t)„4 2 (t)„4 3 (t)) under consideration contain noises in one of the following three forms, (Ci( T a)C2(Tb)Ci( r c)C 2 (Td)), 
(ii(Ta)£,l-i{n)£,i{T c )£,3~i(T d )) and (4( T a)&-i(n)d(Tc)£ 3 _i(Td)), where i= 1,2 and time arguments are dummy vari- 
ables. All 0(<? 3 ) terms in the six moments of the form (*4i(t)*4 2 (t))(*4 3 (t)) under consideration contain the same 
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noises as their corresponding {Ai(t)A2{t)Az,(t)) term. However, in these terms of the form (Ai(t)A2{t))(A?,(t)) 
four noise averages from corresponding terms of the form {A\(t)A2(t)Az(t)) are split into the product of two av- 
erages of two noises. For example, {011^)012^)013 (r)) contains noises in the form (^i(T a )^2(n)^i (t c )^2 ( T d)) whilst 
(ai(T)a2(r))(a^"(r)) contains them in the form (£,i(T a )l;2(Tb)} ( T c)£.2 ( T d)) ■ Using the formula 

^i(r a )Un)Urc)Ur d )) = (6^)6(7*)) {^(^^(tvj)) (3.9) 

+<£l(Ta)&(T e )) (a(n)Ur d )) 

it can be shown that noise expressions in moments of the form {A\{t)A2{t)A^{t)) under consideration factorize. In 
particular, they reduce to the noise expression in the six corresponding term of the form (^4i(r)^4 2 (t)) (A3 (r)). It 
follows that cancellation occurs between the 0(g 3 ) terms in corresponding moments of the form (Ai(t)A2{t)Az(t)) 
and (A\(t)A2{t)) (A3,(t)) under consideration as the two terms are identical. Consequently, all six moments under 
consideration are 0(g i ). They are also typically much smaller than the two 0(g 3 ) moments, (ai^a^MAa^r)) 
and (a^(T)aJ(T)Aa^"(r)), as g <C 1 for realistic systems. To be precise, as all moments are complex quantities, the 
magnitudes of (ai(r)a2(T)Aa3(r)) and (a^ (t)^ (r)Aa^ (r)) are much larger than the magnitudes of the other six 
moments. 

The above results are now more closely related to experiments by considering quadrature phase amplitudes X i ^ i (r). 
In particular, calculations are performed to determine the in principle experimentally observable third order quadrature 
phase amplitude moment (M(r)), where (M(r)) = (AJfi.^ (t)AX2 i 2 (tJA^^j (r)), according to quantum mechanics 
and SED. In quantum mechanics, quadrature phase amplitudes are expressed in terms of creation and annihilation 
operators by the equation 

v a-i exp(-zflj) + aj exp(i^) 

■X-iOt = ■ (6. Li)) 



Using Eq. ( 3.10 ) and operator-positive-P variable correspondences (M(t))qm, the value of (M(r)) for quantum 
mechanics can be expressed as 

(M(r)) QM = J([[A a ,(r)r* + A a +(T)e««'). (3.11) 
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Upon expanding the right hand side of Eq. (3.11), the two lowest order terms in g, (Aai(r)Aa2(7")Aa3(r)) and 
(Aa|(r)Aa|(r)Aa^ (r)), usually dominate. When they do 

(M(t)) qm i[coseRe(ai(T)a 2 (T)Aa 3 (T)) - smQlm{ai{T)a 2 {T)Aa s (T))}, (3.12) 
where 9 = 0\ + 62 + O3. However, w hen c os6 = and Im((ai(r)a2(T)Aa3(r))) = 0(r 4 ) or when sin6 = and 



Re((ai(T)a2(T)Aa 3 (T))) = 0(t 4 ) Eq. ( 3.12 ) is not necessarily true. Such situations can be avoided though because 
O and e are controllable parameters. They are ignored in present considerations. When e is real, the 0(g 3 ) term in 
(ai(T)a 2 (T)Aa3(T)} is also real and so 

x> 1 „ , , , ... , e 2 o 3 cos 6 exp(— 37r) r exp(7r) „ cxp(— jt)-, . . 

(M(t))qm c-cosQ ( ai (r)a 2 (r)A« 3 (r)) ~ y - — ^ U. - 2r - -^ JLL ] ■ (3.13) 



Thus, Eq. (3.13) shows that (M(t))qm is cubic in g, within the domain considered, as shown in Fig. pi 



IV. COMPARISON OF QUANTUM MECHANICS AND STOCHASTIC ELECTRODYNAMICS 

This section compares the predictions of quantum mechanics and SED for the intracavity moment (M(r)). SED 
is a semiclassical theory which adds Gaussian white noise to classical electrodynamics. It describes electromagnetic 
field modes by complex field amplitudes (3. For the nondegenerate parametric oscillator, the set of such amplitudes 
{Pi, /3a j /%} evolves via the equations 
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^ = -7iA + 5/3 2 */3 3 + V7T6 (4-1) 

^ = ~72/3 2 + 3/3^3 + V726 

TT^ = -73/?3 - 3/3i/?2 + V73&, 
or 

where the same time variable as in the quantum case is used and the £'s are independent complex Gaussian white 
noises with the following correlations 

te(n) £*(t 2 )> = S^n - 7a), (4.2) 

where i, j = 1,2,3. The field amplitudes 0x, /3 2 and /3 3 initially have Gaussian fluctuations in their real and imaginary 
parts of variance 1/4. The only nonzero correlations present in these fluctuations are thus 

(AA(0)A/3*(0)) = i, (4.3) 

where i = 1, 2, 3. Initial conditions are (/?i(0)) = (/3 2 (0)) = and (/3 3 (0)) = e. 

The SED prediction for the intracavity moment (M(r)) is {M(t))sed, which is given by the equation 

(M(t)) S ed = (AI 1A (r)AI 2)(j (T)AI 39! (r)>, (4.4) 

where X^g^r) = ( y (3i{r)e~ iei + /3*(r)e iei )/2. It is calculated using a similar iterative method to the one in Section y, 
except that noise terms are now treated exactly instead of iteratively. Zeroth order approximations for this iterative 
method are thus 

Pi 0) (r) = A(0) exp(- 7l r) + / d n expfafo - r)]^i (4.5) 



1"i = 



4°\r) =/3 3 (0)exp(- 7 3T), 



where z = 1,2. Higher order (m + l)*' 1 order approximations are 

4 m+1) W = $ 0) (r) + / n_T dr, exp[ 7i (r, - r)]g^ { n )tff> (r«) (4.6) 



Ti=0 
T3 = T 



f3t +1) (r) = 4°\r) - / dT 3 exp[ 73 (T3 - r)]g(3[ m) (r 3 )^ (r 3 ), 

Jt 3 =0 

where « = 1, 2. 

The lowest order nonzero term in g of (M(t))sed is now found using the same method as for the lowest order 
nonzero term of (M(t))qm ■ Consider the moments of the form (ASi(t)A,B 2 (t)A,B3(t)), where B n (r) is either /3„(r) 
or /3* (t) . The stochastic diagrams required to determine the order of the lowest order nonzero terms of these moments 
are shown in Fig. || Note that noise terms are now 0(g°), instead of 0(g?) as for quantum mechanics. Using the 
stochastic diagrams in Fig. |] it is found that, when 71 = 72 = 73 = 7 and g, Tf <C 1, 

(AA(r)Aft(r)A^(r)) ~ ^-[1 - exp(-3 7 r)] ~ (A#(r)A/$(r)Aft(r))*. (4.7) 

The other six moments of the form (AB 1 AB 2 AB 3 ) are all 0{g 2 ). Thus, (A/3i(t)A/3 2 (7-)A/3|(t)) and 
(A/3*(t)A/3|(t)A/33(t)) dominate these other six moments when g 1 and hence 

I n cos <t> 

(M(t)) S ed -cos$(A/3 1 (r)A/3 2 (r)A/3*(r)) ~ 9 — [1 - exp(l - 7 r)], (4.8) 

4 47 

where $ = 61 + (9 2 — 6*3, when $ 7^ 0. Eq. (h8) shows that (M(t)) sed is linear in g, as shown in Fig. [| This is in 
contrast to the cubic behaviour of (M(t))qm ■ Thus, quantum mechanics and SED predict greatly different values for 
(M(t)) when g -C 1. 

Consideration now is given to the effect of damping strength on the size of the difference between (M(t))sed and 
(M{t))qm ■ Fig. H (a) shows (M)qm and (M)sed as functions of 7 for g = 0.1, r = 1 and e = 1. It indicates 
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that the difference between them is somewhat sensitive to 7, decreasing exponentially with increasing 7 and quickly 
approaching zero. However, for the shorter time r = 0.1, Fig. ||(b) shows that this difference is not as sensitive to 
damping, ft, approximately, only decreases linearly with increasing 7. 

SE P a nd the pos itive- P representation treat fluctuations very differently, as is evident by comparing noise terms in 
Eqs (2.3) and (4.1). This difference in treatment underlies the differences between the two theories' results. Firstly, 

noise terms in the positive-P representation are nonlinear and are scaled by either yjgai or <J ga^ , whilst those in 
SED are linear and are scaled by ^/t7< Secondly, noise terms possess different correlations in the two cases. Thirdly, 
in quantum mechanics no energy fluctuations occur in the vacuum state, whilst in SED {Pi} fluctuates, as does the 
total energy. Assuming quantum mechanics is true, in SED fluctuations in the vacuum lead to an overestimate of 
(M(t)) for small g. 



V. NUMERICAL RESULTS 



The analytic results for (M(t))sed and (M(t))qm in Sections III and IV only include lowest order nonzero terms. 
This leaves the sums of all higher order terms as neglected and these may be significant. For this reason, the validity 
of the analytic approximations are checked by comparison with highly accurate numerical simulation results. 

Numerical simulation methods for stochastic differential equations (SDE's) are both somewhat complex an d not 
widely known. Thus, explanations arc given for the numerical technique used to solve the SDE's in Eqs (2.3) and 
(4.1). Normal ODE techniques such as the Runge-Kutta method cannot be used to solve SDE's as they contain 
discontinuous source terms. Instead, a semi-implicit numerical method |25j| is employed. Only its application to 



Eq. (2.3) is explained as its application Eq. (4.1) is similar. Each of the equations in Eq. (2.3) can be rewritten as 



dxj 
dr 



A(x) + ^B 13 (x)0(r), 



(5.1) 



where x$ is either on or af , for i — 1,2,3, x is a vector whose components are {ai,af}, Ai is the function of x 
formed by the damping and nonlinear terms in the evolution equation for n and bij is a matrix whose elements are 
coefficients of the noise terms {Q} where Q is either £j or £?, for j — 1, 2. The semi-implicit method used determines 

x^ n \ an approximation to x at the midpoint of the interval (r n ,r„ + i). This approximation is found using iteration 
such that the p th order approximation to a component of x~i W W is given by the equation 



» 



(5.2) 



where x\ is the value of Xi at time t„ , At = r n _i — t„, AWj (r„) = rf"^ (r„) At and r„ is the midpoint of the interval 



(j n — i,T n ). The zeroth order approximation to x[^ is given by the equation The approximation to yS n ' 

calculated is then used to generate Ax[ n \ an approximation to the change in Xi over the interval (r„, T re+1 ). This is 
done by solving the equation 



» 



Ax\ n) = A(x (n) )Ar„ 



3 



B y (x (n) )AW,-(f„ 



(5.3) 



Repeated use of Eq. (5.3) determines x+ for successively later and later times and thus solves Eq. (5.1). Two of the 
most important parameters used in the numerical simulations are the step size and the number of stochastic paths 
that are averaged over. The former is always 0.0025 and the latter is O(10 6 ) for most simulations. However, large 
sampling errors necessitated averaging over O(10 7 ) paths for g = 0.1 SED simulations. 

Results from the numerical and analytic simulations of (M(t))sed and (M(t))qm over a range of g and N values, 
where N is the average initial number of pump photons (JV = |e 2 |), are shown in Figs|(j-[| In all cases 9\ = #2 = #3 = 
and relative numerical errors are small. All g = 0.1 analytic results are in agreement with their numerical counterparts. 
However, g = 1 analytic results for N=l and N=10 are not. This disagreement is explained by noting that the analytic 
results are only necessarily valid when g -C 1. 

A number of qualitative trends can be seen in Figs §§ In Figs I and @ (N=l and N=10) the results of SED and 
quantum mechanics are so distinct that they have different signs, with those of quantum mechanics being negative 
and those of SED being positive. This trend only holds for short times (r < 0.07) in Figs @(a) and (b) (N=100). For 
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longer times, SED and quantum mechanics predict the same sign. This trait is consistent with the fact that Figs ra(a) 
and (b) show results for the largest number of photons in the pump mode. SED and quantum mechanics are at their 
most classical level for this case and thus might be expected to differ the least. For constant N Figs [||| also show that 
as g is decreased the results of quantum mechanics and SED become more similar. This occurs because lower g values 
are associated with larger damping to nonlinear coupling ratios and therefore move SED and quantum mechanics 
closer to the classical domain. 



VI. EXTERNAL MOMENTS 



Thus far, only intracavity fields have been considered. However, it is the external fields that leak out of a cavity that 
are observed. In realistic systems, intracavity photons are transmitted through imperfect mirrors into the external 
environment where they are detected. Thus, an external field analogue of (M(t)), (M^ e \t s , 17)), where r s and r/ 
are initial and final measurement times, is calculated according to quantum mechanics and SED to consider what is 
actually observed in the laboratory. 

The first step in calculating the external moment [M^ e \t))qm for quantum mechanics is defining the external 
quadrature phase amplitudes constituting it. This is done within the the context of homodyne detection as quadrature 
phase amplitudes are commonly measured using it. A schematic diagram for balanced homodyne detection is shown 
in Fig. ||. An external signal field flux $i out, where i = 1,2,3, and a local oscillator field flux E{ are incident on a 
50-50 beam splitter BS. An external local oscillator phase variable is represented by Qi. The two field fluxes combine 
and are detected by two photodiodes D + i and D-{. The detected photocurrents are then converted to amplified 

electrical currents whose difference is found. An external quadrature phase amplitude for quantum mechanics X: J 
is defined as this difference yielding, when E t is real, 

^(|) (r) = eA^E^ouT^e^ + $j ot/r ( T )e ifl ~^ (gJ) 

where e is the magnitude of the charge of an electron, A is an amplification factor and r\i is a detector efficiency factor 
for both detectors associated with external field modes denoted by i. In realistic experiments detection occurs over a 
finite period of time and thus 

T=T 'pl$(r) (6-2) 

corresponds to what is observed. Only the 8i = case is considered. Thus, an external moment analogue of the 
intracavity moment (M(t))qm can be defined as 

(M^(r s ,r f )) QM = <r- 3 ;Q r^'dnAxW („)) (6.3) 

l=1 Jr i= r s 

r-Tl=Tf fT 2 =Tf pT 3 =Tf 3 

J Tl =T s Jt 2 =t s Jt 3 =t s i= j 

To calculate (M^ (r s , r/)) qm, the relation between the unknown external output fields that define it and known 
intracavity fields needs to be ascertained. Gardiner and Collett [^(| have formulated an input-output theory which 
relates the two via the equation 



$i out(t) = V2ra,(r) + m (r), (6.4) 

where <!>i in(t) is the input field flux associated with intracavity mode i. All input fields are assumed to be in vacuum 
states. This allows the use of Eq. (5.3) from p6[ , which can be expressed as, in this paper's notation, 

our( Tn )^ i ouT{i~' n+1 ).--&i ouT{T' m )) (6.5) 
= (2rr/2 < f[aj( Tl )aJ( T2 )...aj( T „)]T[a i (r' x).-^)]), 



where T and T are time anti-ordering and time ordering operators respectively. Using Eq. (6.1), the integrand of 



Eq. ( |6.3| ) can be expressed in terms of QiouT and ^\ OUT - It can then be expressed in terms of particular eij(Tj) and 







aj(rj) averages using Eq. ( |6.5| ). In turn, these averages are equivalent to the eight positive-P averages of the form 
(AAi(ti)AA2(t2)AA3(t 3 )}, where Ai is either a, or af. As was determined in Section [II, two of these averages, 
(Aai(n)Aa2(r2)Aa3(r3)) and (Aa^(ri)Ao! 2 f (r 2 )Aag (t 3 )), are of lower order in g than the others and hence dominate 
when g <C 1. Thus, the external field moment of quantum mechanics (M^ e \t s , t/))qm can be expressed as, when 
= 0, where i = 1,2, 3, 



V2(eA V E) 3 r~ 3 / 2 



Tf pT-2—Tf pT^—Tf 



(M^^r^gM ^^ 1 / / / (6.6) 



7"3 = 



(Aai(Ti)Aa2(r 2 )AQ;3(r3)) + (Aa+(r 1 )Aa 2 + (r 2 )Aa+(r 3 )), 

where r\ — r/i and E — Ei, where i = 1,2,3. To simplify the algebra only the r s = case is investigated, so that 
only the moment (M^ (t/))qm(= (^^ E \^t t }))qm) is considered. External fields are only considered for small 
times (rf << 1), and so, to a given order in g, (t"/))qm's lowest nonzero order term in r/ dominates. Hence 

(M' e '(t/))qm can be approximated by its lowest nonzero order term in both g and Tf. Thus, 

(M^(T f )) QM ~ -^^( eA? £) 3 r- 3 /Vf (6.7) 
The SED external moment (M^ (Tf)) sed is now calculated. It is given by the same expression as (M(t/))qM) the 

I — k * ( w\ 

right hand side of Eq. ( p.3| ) (when r s = 0), except that the quadrature phase amplitude operator X^ s '(t), is replaced 
by its SED c-number analogue. This external SED c-number quadrature phase amplitude is defined as, when Ei is 
real, 

y {E) M eAruEiiPi our^e^ + P* OUT (r)e iS >) ( . 

where Piour^) is the output field flux associated with the intracavity field denoted by i. In analogy with Eq. (|6.4|) , 
it is assumed that the SED input-output relation is 



/Wt(t) - V2Tf3 l (T)+p tIN (T), (6.9) 

where Phn(t) is the input field flux for the intracavity mode i. When all input fields are in vacuum states, as is the 
case, Phn(t) is a Gaussian white noise with a self correlation characterized by 

</Wri)#™W)> = 6{Tl 2 T T[) - (6 ' 10) 



A ca lcul ation analogous to the quantum mechanical one earlier in this section can be performed using Eqs (6.S) 
and (6.9) to obtain an expression for (M^ E > (Tf)} sed in terms of particular intracavity averages. When lowest order 
nonzero approximations to these averages are considered, the following result is obtained when 9i for i — 1,2,3, and 
9,Tf < 1, 

(M^(T f )) SED ^ ~ gTf(eAr,E) 3 T- 3 / 2 . (6.11) 

Upon c omp aring Eq. J6.ll ) to result of quantum mechanics in Eq. (3.7), it is seen that the leading order term in g 
in Eq. (6/7) is 0(g 3 ) whilst in Eq. ( 6.11 ) it is O(g). Hence, as was the case for the intracavity moment, quantum 
mechanics and SED predict significantly different results for the observable external field moment (M^ E \Tf)). 



VII. SIGNAL TO NOISE RATIO 



In actual experiments, only finite samples of results are obtained, as opposed to infinite ones. Hence, in practice the 
population means considered thus far are estimated from sample means. These sample means fluctuate from sample 
to sample and thus have signal to noise ratios, which are now determined for small times (r/ << 1). This paper 
focuses on differences between quantum mechanics and SED. Thus, a calculation is performed of the signal to noise 
ratio of the difference between the two theories' external sample moments. First, the noise of the external sample 
moment in quantum mechanics is determined. It is then assumed that the noise of the external sample moment of 
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SED is the same. Noise results are combined with the external moment results of Section VI to produce S(t/), the 
signal to noise ratio of the difference between the two theories' external sample moments. This quantity S(rf) is given 
by 

s( ) = (\(m {E Hrf))sED-(m( E Hr f )) QM \) 



VW=T\(MW( Tf ))sED - (M (E} (r f )) 



QM\ 



V2<j(M( E )( Tf )) 

where n is the number of observations in the sample considered, (m^ E \Tf)) sed and (m' E '(r/))QM are sample averages 
of (M( e \tj)) according to SED and quantum mechanics respectively, and s 2 (A) denotes the sample variance of A. 
The only significant unknown quantity on the right hand side of Eq. (7.1) is a (M( B )(r/)), which is now determined. 
Expressing a(M^) explicitly yields 



v(MW(T f )) = ^(M (E Hrf) 2 )QM-(M^(r f )) 2 QM - (7.2) 

The moment {M^ E \rf))QM was determined in Section[v| and so {M^ E \Tf) 2 )QM is now calculated. In the calculation 
that follows only the Oi = 0, where i = 1,2, 3, and g, r/ <C 1 cases are considered. 

The moment (M^ (t/) 2 )qm can be expressed in terms of external quadrature phase operators as 



3 m=Tf 



(m (£) (t/) 2 )qm = <[r- 3 n / dnAXl E \ n )f) (7.3) 

t=1 Jr t =0 



pT\=Tf r T \— T f PT2—Tf r-r'^—Tj pT^—Tf />T.^—Tf 3 

= r- 6 / / / / / /■ dT 1 dT{dT 2 dT<dT 3 dTl(]lAXl E \T l )AXl E \Tl)), 

Jn^O JtJ^O Jt 2 =0 Jr^=0 Jt 3 =0 Jr^=0 



where X\ '(ji) = X)g =Q ( T i)- The integrand of Eq. ( [7.3| ), which is denoted by K, can be expressed as 

K = f[(AX {E \n)AX^ E \r^ + f( n , t<), (7.4) 

i=l 

where /(t^, t-) is a function which includes terms resulting from coupling between modes. These coupling terms vanish 
when g — and thus are at least 0(g). It follows that K can be be re-expressed as 

K = I]<AXf \ Ti )AX^(ri)) + 0(g). (7.5) 

i=l 

( E) * ( E) 

The moment {AX^ ' {n) AX\ a > {t 1 ^)) is now calculated using a normally ordered approach that has been previously 
employed to solve similar problems p7|, p3|. This method expresses (AX^ \ri)AX^ E (r 8 ')) in terms of normally 
ordered photocurrent averages and then determines these averages. It first defines X\ (r a ), where r a is any r variable, 
as the difference between the amplified electrical currents, XW(t ) and X_ i (r a ), produced by the p/io£ocurrents 
detected at the detectors £>+ i andD_i in Fig. || in Section VI. Using this definition (X- E \T a ) = X^ (r a ) — X_^ (T a )) , 
(AXi E) (Ti)Axl E) (Tl)) can be expressed as 

{Ax\ E \n)Axl E \ri)) = {{Axf}{n) - AX^(n))(^X^(r:) A^(r!))). (7.6) 

Upon expansion, the right hand side of Eq. ( |7.6| ) contains two types of terms, those of the form (X^ (r a )) , where C 

is either + or -, and those of the form (X^, . (Ti)X D ; (t-)), where D is cither + or -. Terms of the form (X^,, (t )) 
are given by the equation 



(X (E }(r a )) = I ^l G g) ( Sl )J°(r a - Sl ), (7-7) 
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where G Ci (si) is a first order Glauber correlation function and J°(r a — s±) is an electrical current pulse produced by 
a single photodetection event. In following previous work P7fl, square electrical current pulses of the form 



J° (a - b) = i 



_ / Aer/T d b<a<b + r d 



(7.1 



are considered in the limit of — > 0, which is taken at some appropriate later stage of the calculation. The Glauber 
correlation function G^\(si) can be expressed as a power series in g and S\ and thus as Ylm'n^Q c mnS"' s r Due to the 
form of J°(r a — si), when r a << 1, as is being assumed, only photodetection events at small times si, contribute to 

" ( E) 

(Xq I {To)) ■ This fact, coupled with the knowledge that only the g « 1 case is considered, means that the n = m = 
term in the power series for Gq Jsi) dominates when J°(r a — si) is nonzero. Hence, upon calculating this dominant 
term by expressing and $J in terms of intracavity field operators, in the limit of large local oscillator amplitude, 



VCi E 2 



where ijc % is a detector efficiency factor for the photodetector Dc i ■ It follows that 



(7.9) 



(7.10) 



" (E) " f_E) I ■ 

Terms of the form pQf/ (t;)^/^)) in Eq. (|7J) can be expressed as 



^Gg)( Sl )jW(r,- Sl )j(°)(r/- Sl ) 



(7.11) 



efelcLs2 (2 



where G c Di (s\,S2) is a second order Glauber correlation function and 5c, d is one when C and D are the same and 
zero otherwise. In the limit of — > 0, 



(7.12) 



to leading nonzero order in g, r» and t 4 '. It is of equal order in g and lower order in Tj and r z ' than the second term 
in Eq. (7.11) and hence is much larger than this second term when it is nonzero as the t^t/ << 1 case is being 



considered. Thus 



(E), ,x. _ ScpjAeEjfricrfin - r/)r 



(7.13) 



From Eqs ( [7.10 ) and ( 7.1 3| ) it can be seen that the single integral terms in (X^(tj)X^(t|)) and 
(X_ i (tj) X_ i (t- )) are of the same order in g and lower order in n and t[ than any other terms contributing to 
(J*Q '(rj)JQ' (r-)) and hence dominate. It follows that 



{xf\n)Xr\Tl))-(AeEf m 5{n 



r/)T, 



(7.14) 



where r/i = ryc^ = and E = E.- L . where i = 1,2,3. As right hand side of Eq. (7.14) is O(g ), 

nf=i(Alf E) (rjJAlf' (t<)) is also O(g ) and hence from Eq. Q, 



if -((Ae^r^n^-r?). 



(7.15) 



Substituting this approximation for K into Eq. (|7.3j) yields 
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3 / 

(M^r/) 2 ) ~ T^((eAE) 2 V Tff[ [ P'^ dndr^n - r[) (7.16) 

i=1 Jr 1= Jr{ 



(eAE) 2 riT f 



where r\ = where i = 1,2, 3. Thus 



*(M^Hr f ))^[^p^f 2 . (7.17) 



Hence, the signal to noise ratio of the difference between the external sample moments of quantum mechanics and 
SED is 

S(r f ) ~ L f . (7.18) 



VIII. REALISTIC SYSTEMS 

Realistic parameter values are now considered to determine if the theoretical difference between SED and quantum 
mechanics could be observed experimentally. In particular, the signal to noise ratio of the difference between the 
sample moments of quantum mechanics and SED 5* is calculated using realistic parameter values for nondegenerate 
parametric oscillators containing the commonly used crystals, silver gallium sclinidc (AgGaSe2) and potassium titanyl 
phosphate (KTP). The non-linear interaction strength G for parametric down conversion is given by (29) 



G " 4// V e V V (8 - 1} 

where V is the cavity volume, 1 is the crystal length and L the cavity length. Cavity and crystal length values of 
10cm are chosen. The cavity volume V is given by the formula V — ir£l 2 L, where Q is the spot size. This volume 
is minimized in order to maximize G and thus the external difference between quantum mechanics and SED. It is 
assumed that the damping constant used to scale time T equals the unsealed damping constant for each mode I\ 
(r = Ti). This common damping constant V is calculated from the formula r = Tx where c is the speed of light 
and T is a mirror transmission coefficient. A T value of T = 0.01 is used. Using the above information, Table [| shows 
realistic parameter values for d e ff, V, G, pump, signal and idler wavelengths, and resulting g and T values. Results 
for {M {E) ) QM and (M^)sed are obtained using Eqs |[7|) and fl6.1l| ) for when rj = 1, r/ = 0.1, E = lO^ 1 / 2 , 
A = 1/e and e = 10 3 . These are displayed in Table 0, which shows that the external results of quantum mechanics 
and SED differ greatly. Due to local oscillator amplification, they are also macroscopically distinct with respect to 
photon number, even though the initial number of intracavity photons is small on average. Another appealing feature 
of the difference between the two theories is that detector efficiencies approach one as photodiodes as opposed to 
photomultipliers are used for detection. Thus, no fair sampling assumptions need to be made. 

The question remains of whether or not the population difference between SED and quantum mechanics could be 



reliably observed in a finite sample of results. To answer it, S is now considered. Figs |10| and 11 show graphs of S 
versus sample size n for AgGaSe2 and KTP for the same parameter values as used in the last paragraph. These show 
reasonable S values and indicate that large sample sizes must be obtained to produce a signal to noise ratio of one, 
the smallest signal to noise ratio required to clearly observe the signal. In particular, sample sizes of 1.8 x 10 13 (KTP) 
and 3.7 x 10 11 (AgGaSe2) need to be obtained to generate a signal to noise ratio of one. An individual observation 
takes a time of the order t = Tf/T = 6.7 x 10~ 9 s and so, assuming minimal time delay between measurements, 
1.8 x 10 13 observations would take about 33 hours and 3.7 x 10 11 observations about 41 minutes. It is conceivable 
that measurements could be taken over both times. Furthermore, as the signal to noise ratio scales as i, higher g 
materials would enable the difference to be observed even more readily. 

IX. DISCUSSION 

It has been shown that there exist a significant, potentially experimentally observable, difference between quantum 
mechanics and SED. Due to local oscillator amplification, this difference can involve macroscopically distinct external 
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fields for the two theories. Thus, it can be considered macroscopic if it is legitimate to include the local oscillators 
as part of the system and not as external measuring apparatuses. The difference is also potentially experimentally 
observable, as a realistic system and state are considered and is present at realistic parameter values. The system 
is practical as parametric oscillators and balanced homodyne detection are widely used, and damping is included. 
The state is realistic as the initial intracavity coherent state can be approximated well by a laser. It follows that the 
difference can be seen as providing the basis for an experimentally achievable macroscopic test of quantum mechanics 
against one local hidden variable theory (SED). Such a test is significant as all experimental tests of quantum mechanics 
against local hidden variable theories to date have been microscopic. It is true that many macroscopic tests have been 
proposed, but most of them consider highly idealized states or systems that are not currently able to be experimentally 
implemented. In particular, many of them do not consider damping, even though it is known to rapidly destroy the 
correlations of quantum mechanics present in Schroedinger cat p2[ and other entangled states. The calculations 
in this paper do include damping and show that the difference between SED and quantum mechanics is not overly 
sensitive to it. Most importantly, it remains for realistic damping values. The test proposed in this paper can be seen 
as being in the novel and largely unexplored domain of macroscopic experimental tests of quantum mechanics. 

Even if the local oscillators are not included as part of the system investigated, the external difference between 
quantum mechanics and SED is still at least mesoscopic as average initial pump photon numbers up to 10 6 are 
considered. From this perspective, the difference is still distinct from many earlier microscopic ones known to exist 
between quantum mechanics and all local hidden variable theories. It is also, perhaps, more surprising than some of 
them as it occurs in a larger particle number system. 

Two noteworthy features of the external difference between quantum mechanics and SED are that it involves contin- 
uous variables and high efficiency detection. That it involves continuous variables is significant because most previous 
differences between quantum mechanics and local hidden variable theories have involved discrete ones. Furthermore, 
it is, perhaps, more surprising that a difference between quantum mechanics and a local hidden variable theory can 
be found for continuous variables as continuous variables are more closely related to classical ones (which are all 
continuous) than discrete ones. Low detector efficiency forms the basis of a significant loophole in most tests between 
quantum mechanics and SED to date (3^] . The use of photodiodes for detection in the scheme discussed means that 
such a loophole is avoided. 

The calculations in Section VIII show it is difficult to observe the external difference between quantum mechanics 
and SED. This is mainly a result of small experimental nonlinearities. They cause few signal and idler photons to be 
created and thus the experimental signal is weak relative to its noise. For small enough measurement samples, SED 
results cannot be clearly distinguished from those of quantum mechanics. This fact is consistent with the knowledge 
that SED reproduces many features of quantum mechanics. However, it is a distinct theory and does differ from 
quantum mechanics in particular cases, as this paper has shown. 

The external difference between quantum mechanics and SED would be easier to observe if larger nonlinear cou- 
pling constants were used. These could be achieved by using organic nonlinear crystals such as iV-(4-nitrophcnyl)-L- 
prolinol (NPP) Q]. However, phase matching would be difficult with such crystals. In addition, they are typically 
only transparent within a small frequency range. Alternatively, higher nonlinearities could be achieved by using 
Josephson-parametric amplifiers ]35|) , which can have even larger nonlinearities than organic nonlinear crystals. An- 
other possibility, in the area of atom optics, is to to utilize BEC nonlinear effects, in which atom-molecule coupling is 
induced through photon-associaton ||(| . 

To conclude, this paper compared particular moments of quantum mechanics to those of SED for the nondegenerate 
parametric oscillator. Both internal and external moments were considered and an analytic iterative technique showed 
them both to be cubic in the system's nonlinear coupling constant for quantum mechanics and linear for SED. 
Numerical simulations were performed to check the approximate intracavity analytic result and were in agreement 
with them when the system's nonlinear coupling constant was much less than one. Realistic parameter values were 
considered and it was shown that the external sample difference between SED and quantum mechanics had a small 
signal to noise ratio in typical parametric oscillators. The presence of intense local oscillators means that the results 
could be seen as providing the basis for a macroscopic experimental test of quantum mechanics against SED. 
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Figure 1 : Pope et al. 




FIG. 1. The three basic classes of stochastic diagrams, (a) initital value term, (b) noise term and (c) nonlinear term. A 
stochastic diagram (d) representing a higher order term. 
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Figure 2:Pope et al. 




FIG. 2. Stochastic diagrams representing the lowest order nonzero terms required to determine the intracavity moment of 
quantum mechanics (M(t))qm for (a) ai(r),a\(r),a2(r) anda^r) and (b) ck3(t), c4(t), (0:3 (r)) and (al(r)). 



Figure 3: Pope et. al. 
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FIG. 3. Analytic results for (M(t))qm (solid line) and {M(t))sed (dotted line) versus scaled time r for N=l, g=l, 7 = 1 
and cos = cos $ = 1. 
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Figure 4: Pope et al. 




FIG. 4. Stochastic diagrams representing the lowest order nonzero terms required to determine the intracavity moment of 
SED (M(t)) S ed for ft, where i = 1, 2, 3. 
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Figure 5a: Pope et. al. 
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Figure 5b: Pope et. al. 
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FIG. 5. Results for {M(j))qm (dotted lines) and {M(y))sED (solid lines) as a function of the damping constant 7 for 
g = 0.1,7V = 1, cosB = cos* = and (a) r = 1, (b) r = 0.1. 



Figure 6a: Pope el. al. 
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Figure 6b: Pope et. ; 



FIG. 6. Numerical and analytic results for {M(t))qm and {M(t))sed for N = 1, 7 = 1 and cosG = cos$ = 1. (a) Analytic 
results for {M(t))sed for g = 0.1 (solid line) and g = 1 (dotted line), and numerical results for (M(t))sed represented by 
dots with associated error bars for i) g = 1 and ii) g — 0.1. (b) Analytic results for {M(t))qm for g = 0.1 (solid line) and g = 1 
(dotted line), and numerical results for {M(t))qm represented by dots with associated error bars for i) g = 1 and ii) g — 0.1. 



Figure 7: Pope el. al. 




FIG. 7. Numerical and analytic results for {M{t))qm and {M(t))sed for N = 10, 7 = 1 and cos0 = cos$ = 1. Analytic 
results for (M(t))sed are indicated by solid lines for i) g = 1 and iii) g = 0.1. Numerical results for {M(t))sed are indicated 
by dots and with associated error bars for ii) g = 1 and iii) g = 0.1. Analytic results for {M(t))qm are indicated by dotted 
lines for v) g = 1 and iv) (/ = 0.1. Numerical results for {M(t))qm are indicated by dots with associated error bars for v) g = 1 
and iv) (? = 0.1. 
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Figure 8a: Pope et. al. 
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Figure 8b: Pope et. al. 
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FIG. 8. Numerical and analytic results for (M(t))qm and {M{t))sed for N = 100, 7 = 1 and cos 6 = 1 = cos$ = 1. (a) 
Numerical results for {M(t))sed represented by dots with associated errorbars for ii) g — 0.1 and iii) g = 1, and analytic 
results for {M(t))sed for g = 0.1 (dotted line) and g = 1 (solid line), b) Numerical results for (M(t))qm represented by dots 
with associated errorbars for i) g = 0.1 and ii) g = 1, and analytic results for (M(t))qm for <? = 1 (dotted line) and g = 0.1 
(solid line). 
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Figure 9: Pope et al. 
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FIG. 9. Schematic diagram for balanced homodyne detection. Photodetectors are labeled by D+i and D-i, for i = 1,2,3, 
BS is a beamsplitter, 8i is a local oscillator phase variable, Ei is a local oscillator amplitude and <&i is an external signal field 
operator. 
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FIG. 10. Plot of signal to noise ratio of the difference between the sample moments of quantum mechanics and SED S versus 
sample size n for AgGaSe2 for g=8.3 x 10~ 3 , Tf — 0.1 and r\ = 1. 
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Figure 1 1 : Pope et. al 
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FIG. 11. Plot of signal to noise ratio of the difference between the sample moments of quantum mechanics and SED S versus 
sample size n for KTP for g=1.2 x 10~ 3 , 77 = 0.1 and 77 = 1. 
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CRYSTAL 


d ef f(pmV r ) 


V(m 3 ) 


A 3 ~ A 2 /2 ~ Ai/2(>m) 


G(0 


r^- 1 ) 


9{=$) 


AgGaSe2 


33 (at A = 2.1/im) [ 


|| 1.0 x l(T a 


1.4 


1.3 x 10 b 


1.5 x 10' 


8.3 x W-- 1 


KTP 


7.2 (at A = 2.3^m) 


U 1.6 x 10" 9 


1.6 


1.8 x 10 4 


1.5 x 10 7 


1.2 x 10" 3 



TABLE I. Table of realistic values for the parameters d e //, volume V, pump, signal and idler wavelengths A3, A2 and Ai, 
nonlinear coupling constant G, damping constant V and nonlinear coupling constant to damping ratio g for the nonlinear 
crystals AgGaSe2 and KTP. 
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CRYSTAL 


(M^(r f )) 


{MW(Tf)) S BD 


AgGaSe 2 


-2.9 x 10 8 


1.3 x 10' J 


KTP 


-8.8 x 10 5 


1.8 x 10 8 



TABLE II. Table of the external moments of quantum mechanics and SED for AgGaSe 2 and KTP. 
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